# -*- coding: utf-8 -*-
"""
Created on Wed Jun  1 07:06:54 2024

@author: zhao
"""

from scipy.special import yv # 导入 诺伊曼函数（  第二类  Bessel函数）  # yn 是第二类整数阶Bessel函数
import numpy as np
import matplotlib.pyplot as plt
import scienceplots
x=np.linspace(0.01,10,100)
plt.style.use('science') # 这个需要单独安装
plt.figure(1,dpi=300)  
for n in (-1,-2,-3):    
    plt.plot(x, yv(n,x),linewidth = '0.7', label=r'$\nu $='+str("{:.2f}".format(n))) 

plt.plot(x,np.zeros(100),linewidth = '0.1',color='black') 
plt.title('Neumann function')
plt.xlabel(r'$x$')
plt.ylabel(r'$N_{\nu} (x)$')
plt.ylim(-2,2)
#plt.legend(loc='upper right')
plt.legend()
name='Neumann_负.png'

plt.savefig(name,dpi=300)
plt.show() 